Semiclassical approach to universality in quantum chaotic trans- 
port 
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Abstract -The statistics of quantum transport through chaotic cavities with two leads is en- 
coded in transport moments M m = Tr[(tH) m ], where t is the transmission matrix, which have a 
known universal expression for systems without time-reversal symmetry. We present a semiclas- 
sical derivation of this universality, based on action correlations that exist between sets of long 
scattering trajectories. Our semiclassical formula for M m holds for all values of m and arbitrary 
number of open channels. This is achieved by mapping the problem into two independent com- 
binatorial problems, one involving pairs of set partitions and the other involving factorizations in 
the symmetric group. 
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Introduction. — A remarkable fact about electronic 
transport in ballistic systems with chaotic classical dy- 
namics is that they display universal properties like, e.g., 
conductance fluctuations and weak localization [1] . These 
properties are well described by random matrix theory 
(RMT) m which the system's details are neglected 

and its S matrix is modeled as a random element from an 
appropriate matrix ensemble, specified only by the sym- 
metries present. It has always been a central problem 
to derive such universal results from a semiclassical ap- 
proximation in which classical dynamics is taken into ac- 
count. In this formulation, quantum universality must 
emerge as a result of action correlations that are present 
in any chaotic system: in essence, the existence of long 
scattering trajectories that have nearly the same action 
and interfere constructively in the semiclassical limit [5]. 
The same relation exists between long periodic orbits and 
spectral statistics of closed systems [6]. 

Consider a cavity with two leads attached, with Ni and 
A?2 open channels supported in each lead. This is taken 
as a model of actual situations involving semiconductor 
quantum dots. Let t be the transmission block of the uni- 
tary S matrix. The transport moments M m = Tr[(tH) m ] 
carry much information about scattering through the sys- 
tem. If time- reversal symmetry is broken by a strong mag- 
netic field, the RMT calculation of these quantities can be 
carried out for arbitrary values of Ni and A*2 [7]. However, 
semiclassical derivations that agree exactly with RMT pre- 



dictions are available only for the first two moments [8,9 , 
which are related to the conductance and shot-noise of the 
cavity. Higher moments have been treated perturbatively 
in 1/N, where N = Ni + N2, but only to leading [10] and 
next-to-leading [TT] orders. 

In this article we remedy this situation, providing a 
semiclassical derivation which is valid for all transport mo- 
ments and to all orders in perturbation theory, for broken 
time-reversal symmetry. Starting from a diagrammatic 
formulation, we show how the diagrams that are rele- 
vant for M m , which involve scattering trajectories between 
leads, may be obtained from closed diagrams involving 
only periodic orbits. Then, we map them into certain fac- 
torizations of permutations in the symmetric group. The 
factorizations required are different from the ones in |10j 
(see also [12]) and apparently have not been considered 
before. 

Semiclassical Transport Moments. In the semi- 
classical limit h — > 0, N — > 00, the matrix elements of t 
may be approximated [13] by 
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where the sum is over trajectories starting at incoming 
channel i and ending at outgoing channel o, S-y is the 
action of trajectory 7 and A 7 is an amplitude related to 
its stability. The prcfactor contains the Heisenberg time 
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Tij, which equals N times the classical average dwell time. 
Expanding the trace, transport moments become 
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The sum involves two sets of m trajectories, the 7's and 
the cr's. A 1 — Ylj A lj is a collective stability and 5 7 = 
Y^j is a collective action, and analogously for cr. Most 
importantly, the structure of the trace implies that these 
two sets of trajectories connect the channels in a different 
order, and we can arrange it so that 7^ goes from ij to Oj, 
while <Tj goes from ij to Oj + \. 

The result of ([2]) is in general a strongly fluctuating func- 
tion of the energy, so a local energy average is introduced. 
When this averaging is performed in the stationary phase 
approximation, it selects those sets of cr's that have almost 
the same collective action as the 7's. In the past 10 years 
[14] it has been established that the mechanism behind 
these action correlations is the existence of encounters: 
each cr must follow closely a certain 7 for a period of time, 
and may exchange partners when some of the trajectories 
come together. A region where £ pieces of trajectories run 
nearly parallel and £ partners are exchanged is called an 
^-encounter. The two sets of trajectories are thus nearly 
equal, differing only in the negligible encounter regions. 



In particular, this implies A 1 A* 



This theory has 



been presented with great detail in [T51ITC] . 

We show two examples of correlated sets of trajectories 
in Figure la,b, both contributing to M3. Naturally, this is 
a simplification: in a realistic chaotic system these trajec- 
tories would be long and extremely convoluted. The first 
example, Fig. la, has a triple encounter, within which the 
cr's are represented by dashed lines. Before the encounter 
o~j is indistinguishable from 7^; after the encounter it be- 
comes indistinguishable from 7j+i, i.e. the changing of 
partners inside the encounter happens in such a way that 
CTj starts at ij and ends at Oj+i . Notice that if the en- 
counter was arranged differently, with <jj pairing up with 
7j_i, this would not lead to an acceptable contribution 
to J2]). The second example, Fig. lb, has two double en- 
counters. In this case each a ends up along the same 7 it 
started with. It is only acceptable as a contribution to ^ 
when all outgoing channels coincide. 

The present work will apply exclusively to systems for 
which the dynamics is not invariant under time-reversal. 
This is because we shall not consider situations where a a 
trajectory runs in the opposite sense with respect to a 7 
trajectory. This is a significant restriction on the possible 
correlated sets, and considerably simplifies the treatment. 

A correlated pair of trajectory sets may be represented 
by a diagram with a structure. A diagram is a topolog- 
ical entity, a graph where encounters become vertices of 
even valence and the pieces of trajectories leading from 
one encounter to another (or to leads) become edges [TO] . 
(Channels are formally also vertices, but they should not 
be confused with encounters, so we will keep calling them 




Fig. 1: Two (pictorial) correlated sets contributing to M3 are 
shown in a) and b). The diagrams that correspond to them 
are shown in c) and d), respectively. The structure of a) and 
c) is such that the cr trajectories implement the permutation 
7T = (123) on the channel labels. A different structure could 
be chosen so that n — (132), but in that case c) would be 
considered a pre-diagram, and some channels would have to 
coincide to make it a diagram. Coinciding channels appear in 
b) and d), where only one structure is possible. 



channels.) The structure is a prescription for walking in 
the graph, i.e. it specifies how the cr's change partners at 
the vertices. Examples are shown in Figure [TJ:,d. 

It is known [TB] that the contribution of a diagram 
to transport moments consists of (— l) v N C(N±, N 2 ), 
where V is the number of vertices, E is the number of 
edges and C(N 1: N 2 ) is the contribution due to the chan- 
nels. This factorizes into two parts. First, there is the 
number of ways they can be assigned among the possible 
ones existing in the leads. If the number of distinct in- 
coming and outgoing channels are m\ and m 2 , then this 
is 

Nx\ N 2 \ 



B(mx,m2) 
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Second, when there are k different cr's starting and ending 
at the same channels, there is a fc!, since they could be 
labeled in that many different ways. 

The diagram in FigQJ, contributes as —B(3,3)/N 5 to 
M3. On the other hand, the one in FigQ}) contributes, 
also to M 3 , as 2B(2,2)/N 5 , because a 2 and cr 3 start and 
end at the same channels, so there is freedom in choosing 
which is which. 

It was shown in [TU] that to leading order in 1/N all 
transport moments M m are determined by diagrams that 
have the topology of a tree. This simplifying feature al- 
lowed a recursive approach that lead to an explicit solu- 
tion in agreement with RMT for all m. Later, diagrams 
determining the first corrections were constructed [11] by 
grafting trees on base diagrams with more complicated 
topology. We consider a different method, inspired by the 
treatment of Mi and M 2 in [16]. The crucial point is a 
relation between transport diagrams and certain closed di- 
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Fig. 2: a) The pre-diagram corresponding to the diagram 
in FigQJl. The permutation induced by the cr's is the identity, 
7r = (1)(2)(3). b) Introducing an extra vertex, the pre-diagram 
is turned into correlated periodic orbits. A schematic represen- 
tation of the orbits is shown in FigUJ 

agrams which have no channels but have one extra vertex. 

Pre-diagrams. — We now describe our method, 
which draws on several combinatorial ideas. Details will 
be presented in a future publication. Let S m denote the 
symmetric group of order m, i.e. the set of all permuta- 
tions of m symbols. We denote by l m the identity in S m 
and by c m = (12 ■ • • m) the complete cyclic permutation 
with increasing elements. Note that c m is the permutation 
the cr's are required to perform on the channel labels, i.e. 
they should take ij to Oj+x- On the other hand, the 7's 
perform l m by taking ij to oj. 

Let us introduce a class of graphs we call pre-diagrams, 
which have m incoming and m outgoing channels, assumed 
all distinct. We call m the order of the pre-diagram. Tra- 
jectory 7^ is still required to go from ij to Oj, but Gj is 
allowed to go from ij to any outgoing channel. This means 
that in a pre-diagram the permutation performed by the 
cr's on the channel labels is generally different from c rn . 
Let 7T denote this permutation. We present an example 
in Figure [2^,. Clearly, to every given diagram with given 
structure we can associate a unique pre-diagram (obtained 
simply by ignoring coincidences between channels). 

There are therefore two important differences between 
a diagram and its pre-diagram. First, in a diagram trajec- 
tory Oj start at ij and end at Oj+i, while in a pre-diagram 
it can end anywhere and the set a may implement any per- 
mutation tt on the channel labels. Second, there may be 
coinciding channels in a diagram, while they are all differ- 
ent in a pre-diagram. The idea is that true diagrams can 
be obtained from a pre-diagram by means of those coinci- 
dences among channels that make 7r effectively equivalent 
to c m . For example, if tt = l m then several coincidences 




Fig. 3: The partitions in this graph constitute a coincidence 
pair admissible under the permutation n = (12) (34). This is 
because the bipartite graphs induced by connecting the blocks 
according to n or according to (1234) are the same. The mul- 
tiplicity of this graph is 2, because the trajectories represented 
by the double bond admit two possible labelings. 

are needed, such as all incoming or all outgoing channels. 
In order to count all diagrams, we must count all possible 
coincidences for all possible pre-diagrams. 

Admissible pairs of partitions. — Let us formalize 
the above ideas. Coinciding channels induce partitions of 
the set {l,...,m} on the incoming and outgoing leads. 
Suppose two sets of m points arranged vertically side by 
side. Given a partition L of the left set and a partition R 
of the right set, identify the points according to the blocks 
of the partitions. Given a permutation tt € S m , draw lines 
going from the left points to their images under tt on the 
right. If L and R have mi and mi blocks respectively, this 
produces a bipartite graph Q{L, R; tt) with mi vertices on 
the left and mi vertices on the right. An example is shown 
in Figure[3l where the left partition is L = {{1}, {2, 3, 4}}, 
the right one is R = {{1, 2, 4}, {3}} and the permutation 
is tt = (12)(34). 

We say that L and R form a 7r-admissible coincidence 
pair if Q(L,R;tt) = Q(L, R; c m ). Physically, this means 
that the permutation induced by the cr's is 'effectively' 
equal to c m , and the diagram contributes to ([2]). Let 
A4(L, R;tt) be the incidence matrix of G{L,R;tt). This 
means that Mjk is the number of bonds going from block 
j on the left to block fc on the right. The incidence matrix 

of the graph in Fig(3]is ( ^ J j ■ As we have discussed, 

each multiple bond gives rise to a factorial because the 
a trajectories could be interchanged. We thus define the 
multiplicity of the pair (L, R) as 

fx(L,R;7r) = Y[[M(L,R;ir)} jk \ (4) 

Let A(tt ,mi,mi) be the set of all tt- admissible coinci- 
dence pairs with mi blocks on the left and mi blocks on 
the right. The function 

/(7r,mi,TO 2 )= ^2 h{L,R;tt) (5) 

(L,R)eA(TT.m 1 ,m 2 ) 

counts, with the correct multiplicity, diagrams with mi 
distinct channels on the left and mi distinct channels on 
the right, which have pre-diagram with permutation tt. 
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Once it is known, we can compute 

F(n,N 1 ,N 2 )= f(ir,m 1 ,m 2 )B(m 1 ,m 2 ), (6) 

mi ,m2 

where B(m 1 ,m2) is given by The quantity above is 
the total contribution to the sum © associated with pre- 
diagrams for which the permutation performed by the ct's 
on channel labels is it. Notice that since it <E S m this 
depends implicitly on m. 

We therefore meet our first combinatorial problem, 
which is to obtain the function F(w, N\, N 2 ), or at least 
make it computable via a generating function or a re- 
currence relation. This requires understanding the set 
.4.(71", mi, which seems to be a very complicated com- 
binatorial problem. Once we have F(tt, Ni,N 2 ) under con- 
trol, we are left with the problem of finding all possible 
pre-diagrams. 

Relation to correlated periodic orbits. We now 

turn to the question of counting pre-diagrams. First, we 
associate pre-diagrams with certain collections of corre- 
lated periodic orbits. Let a be a single periodic orbit and 
let /3 be a set of periodic orbits that is correlated with a, 
i.e. /3 differs from a only in encounters. We show an ex- 
ample in Figure HI where one periodic orbit is correlated 
with three others, in a situation with two 2-encountcrs and 
one 3-encounter. 

Given a and /3, suppose we 'cut open' an m-encounter. 
This produces 2m endpoints, m of them corresponding to 
the 'beginning' of trajectories (leaving the encounter) and 
another m to 'ending' of trajectories (arriving at the en- 
counter). We interpret them as incoming and outgoing 
channels, respectively. Then, we choose one of the incom- 
ing channels to be i\, and use a to label all channels in 
sequence: the piece of a that starts in ij (and necessarily 
ends in Oj) becomes jj, while the piece of /3 that starts 
in ij becomes Uj . This obviously produces a pre-diagram, 
with a certain permutation n to be determined. See Figure 
2b. 

There are two caveats in the procedure outlined above. 
The first is that some of the /3's might not participate 
in the encounter we chose to cut open, and would not 
become scattering trajectories. This is remedied simply 
by demanding that in our correlated pair a, /3 there must 
an m-encounter in which all /3's participate. Second, even 
though any pre-diagram related to M m can be produced in 
this way, some of them could be produced more than once. 
We show the simplest possible such example in Figure 5: 
the same pre-diagram could be produced having three /3's 
or only one /3. However, the structure of the encounter 
that was opened is necessarily different. In order to avoid 
this overcounting, we have at our disposal the possibility 
to select the permutation experienced by the cr's inside 
the encounter that was opened. We shall make use of this 
right after Eq.flTJ). 

In |15j . pairs of correlated periodic orbits were associ- 
ated with factorizations of permutations. We summarize 




Fig. 4: (color online) Schematic representation of the corre- 
lated periodic orbits in Fig. 2(b), one depicted with solid line 
and three with dashed lines. With this choice of labels the 
associated factorization equation in the symmetric group S7 
is (1234567) = (1)(264)(375) • (125)(36)(47). Removing the 3- 
encounter we produce the pre-diagram in Fig.2(a), such that 
trajectory 71 has no stretches while trajectories 72 and 73 in- 
clude stretches {3,4} and {5,6}, respectively. 

the idea. Suppose some correlated pair a, (3. Label the 
encounter stretches in such a way that the end of stretch j 
is followed by the beginning of stretch j + This produces 
the permutation eg, where E is the number of stretches, 
acting on the 'exit-to-entrance' space (it goes from the 
exit of an encounter to the entrance of another one). A 
variant of this construction is shown in FigUl where a is 
represented by a solid line. 

The orbits behave differently at the encounters (the 
'entrance-to-exit' space). At any encounter a corresponds 
to the identity permutation, since it takes the entrance of 
a stretch to the exit of the same stretch. On the other 
hand, /3 may be represented by a non-trivial permutation 
P, whose number of cycles is equal to V, the number of 
encounters. In the example shown in Figure 4 we have 
P = (125) (36) (47). The product c E P = Q, acts on 'exit- 
to-exit' space, leading from the exit of an encounter to the 
exit of another one. It must be a single cycle if there is 
a single /3. Since ce is fixed, the total number of corre- 
lated pairs equals the number of solutions in S E to the 
factorization ce — QP -1 . 

In our problem we have a single periodic orbit a which 
is correlated to a set of any number of periodic orbits, /3. 
It is possible to adapt this mapping into permutations to 
aid us in the enumeration of pre-diagrams. We have to 
remove the condition that the permutation Q be single- 
cycle, to allow for more than one /3. In the example shown 
in Figure 4 there are three different /3's, and we have Q = 
(1)(264)(375). We must also ensure that all /3's take part 
in the encounter we open, which we convention to be the 
first one. Finally, we have to make sure that when we 
produce the pre-diagram, the permutation implemented 
by a on the channel labels is given by tt. When we remove 
the encounter involving stretches {1, 2, 5} in Figure 4, we 
produce a pre-diagram for which tt is the identity in S3. 
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Fig. 5: (color online) a) Schematic representation of a peri- 
odic orbit (solid line) which is correlated with three periodic 
orbits (dashed lines). Removing the encounter we produce a 
pre-diagram with permutation n — (1)(2)(3). Notice that the 
permutation inside the encounter, (132) = (321), is decreas- 
ing, b) Schematic representation of a pair of correlated peri- 
odic orbits. Removing the encounter would produce the same 
pre-diagram as a), but we discard this situation because the 
permutation inside the encounter, (123), is not decreasing. 



Some factorizations of permutations. Let {P} 

denote the set of integers which are not fixed points of the 
permutation P. Let Pi denote the first cycle of P, the one 
that contains the element '1'. Given a set s, we define P| , 
the restriction of P to s, to be the permutation obtained 
by simply erasing from the cycle representation of P all 
symbols not in s. For example, (123)1^ 3 \ = (13). 

We also define a slightly more involved operation we 
call reduction. Given a permutation Q and a set s, 
the reduction 1Z S [Q] is obtained by first restricting Q 
to s and then making each element as small as possi- 
ble keeping positivity and relative order. For example, 
^{2,3,4,6} [(264)] = (143). This is found as follows. First, 
because we are restricting to {2,3,4,6}, we must write 
explicitly the fixed point: (264) (3). Then reduction leads 
to (143) (2). Finally, we may omit again the fixed point. 
Another example: 7?.{ 2 ,3,4} [(264)] = (13). 

Suppose we have correlated orbits described by the 
equation ceP = Q- The set of elements involved in the 
first encounter is {Pi}, assumed to have m elements. In 
the example of Fig 0] this is {1,2,5}. The 7 trajectories 
start and end at this encounter and, by construction, visit 
these elements in increasing order, i.e. they implement 
a permutation which is simply Ce\j Pi \- This is (125) in 
FigH 

We must determine what is 7r, the permutation induced 
by tr on those labels. First, we take account the permu- 
tation Q and reduce it to the appropriate space, Q\r Pl y- 
In FigH] this is (1)(2)(5). This acts on exit-to-exit space, 
i.e. it takes incoming channels to incoming channels. We 
must therefore multiply it by the inverse of Pi in order to 
reverse the permutation effected inside the first encounter. 
The result, Q\{ Pl y Pi _1 j takes incoming channels to out- 
going channels. In Figf4]this is also (125), just like for the 
7's. 

At this point, we have the permutations implemented by 
both 7 and a on the channel labels. This first is Ce\s Pi \ 



and the second is Q|r Pl i P x ■ The permutation 7r corre- 
sponds to measuring the second one with respect to the 
first one. In other words, we must carry out a change of 
coordinates of sorts. We therefore multiply both quan- 
tities by the inverse of the first. This turns the second 
c i? 1 |{p 1 }- m our example of FiglU 
The permutation n G S m is finally 
obtained after reduction: 



one into Q\[ Pl \ Pi 
this gives (1)(2)(5). 



n 



{Pi} 



Q\ 



(7) 



As we have mentioned, some pre-diagrams can be pro- 
duced more than once. A simple example is shown in 
Figure [5] In order to avoid this overcounting we may im- 
pose a convention on the permutation experienced by the 
ct's inside the encounter which is cut open. A very con- 
venient choice is to demand that Pj~ = ce\{ Pi }, which 

means that Pj -1 is increasing (i.e. its elements are ordered 
increasingly) or, equivalently, that Pi is decreasing (i.e. 
its elements are ordered decreasingly, as in FigJSK). This 
leads to a much simpler expression for the permutation 7r: 



n = K{Pi} [Q] 



(8) 



In particular, the number of cycles of tt equals the number 
of individual periodic orbits in the set (3. 

All in all, the combinatorial problem that needs to be 
solved is the following. One must find the number of so- 
lutions in Se, let us denote it by S(to, 7t, E, V), to the 
factorization equation ce = QP^ 1 which satisfy several 
conditions: i) P has V + l cycles; ii) P^ 1 is increasing and 
of size m; iii) all cycles of Q have at least one element in 
common with Pi; iv) the reduction of Q to {Pi} is equal 
to a given tt. 

We therefore meet our second combinatorial problem, 
which is to obtain the function S(m, 7r, P, V), or at least 
make it computable via a generating function or a recur- 
rence relation. As we have seen, this requires handling 
a very involved factorization problem in the symmetric 
group. 

Conclusions. — We have shown that the scmiclassical 
calculation for transport moments results in 



M m = J2 F(Tr,N u N 2 )Y,Z(m,Tr,E,V) 



E,V 



N E-\ 



(9) 



where the quantities involved have been defined in the 
text. This expression is exact, i.e. valid for arbitrary num- 
bers of channels. It is easy to implement it in symbolic 
packages and verify that it indeed reproduces RMT re- 
sults, as far as it can be checked. Unfortunately, at present 
both combinatorial problems presented here remain open, 
so that a closed form expression for M m is beyond reach. 
However, we stress that the semiclassical approach is not 
intended as a computational tool. Its merit is in revealing 
the dynamical origin of universality. 
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Exact semiclassical results are somewhat surprising in 
the regime of finite number of channels, specially for high 
moments, because the number of channels in a lead is small 
when its width is of the order of the electronic wavelength, 
when diffraction effects are expected to be important [T7] . 
It is not clear what role is played by diffraction in bring- 
ing about universal statistics in chaotic geometries, and 
how it would affect the semiclassical description presented 
here. Another point that remains to be addressed are the 
corrections arising when the system's Ehrcnfcst time Te 
is larger than the typical dwell time To in the cavity (we 
have assumed Te *C Td , which is when RMT universality 
is expected). Some progress has been made in that direc- 
tion [18], but only to leading order in 1/N. These issues 
deserve further study. 

In summary, we have presented a semiclassical deriva- 
tion of the universal statistics of quantum transport in 
chaotic systems. This was achieved by reducing the prob- 
lem to two independent combinatorial questions. One of 
them involves certain pairs of set partitions and is em- 
bodied in the function F(ir, N\, N 2 ). The other involves 
certain factorizations in Se and leads to the function 
S(m, 7r, E, V). We hope that the present article will draw 
attention to the profound relationship between quantum 
chaos and combinatorics. Similar ideas can be developed 
for closed systems, to provide all spectral correlation func- 
tions. Work in this direction is underway |19j . 

* * * 

This work was supported by FAPESP and CNPq. After 
it was completed I learned that Berkolaiko and Kuipers 
also found a semiclassical derivation of M m valid for all m 
and all channel numbers [20] . 
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